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Abstract 

We study the possibly existing anisotropy in the accelerating expansion Universe with the 
Union2 Type la supernovae data and Gamma-ray burst data. We construct a direction-dependent 
dark energy model and constrain the anisotropy direction and strength of modulation. We 
find that the maximum anisotropic deviation direction is (/, b) = (126°, 13°) (or equivalently 
{I, b) = (306°, —13°)), and the anisotropy level is go = 0.030;|]:[]!^[j (obtained using Union2 data, 
at Icr confidence level). Our results do not show strong evidence for the anisotropic dark energy 
model. We also discuss potential methods that may distinguish the peculiar velocity field from the 
anisotropic dark energy model. 
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1. INTRODUCTION 



The cosmological principle has played an important role in modern cosmology [T]. It 
tells us that our universe is homogeneous and isotropic on large cosmic scale, which is 
consistent with currently observational data sets such as the cosmic microwave background 
(CMB) radiation data from the Wilkinson Microwave Anisotropy Probe (WMAP) [2-0]. Up 
to now, current astronomical observations are still in good agreement with ACDM model 
generally |5]. 

Despite the fact that the concordance cosmological model [ACDM model) is confirmed 
by many observational data, the model is also challenged by some observations [HI E] (see [S] 
and references therein for more details). So, it is necessary to revisit the assumption of 
cosmological principle. As more Type la supernovae data [9l [10] and high-redshift Gamma- 
ray burst (GRB) data are released [HI [12], one can use the supernovae data and the GRB 
data to detect the anisotropic direction of cosmic expansion. 

Theoretically, a variety of cosmological models may give rise to cosmic anisotropy. For 
example, peculiar velocities may generate dipole-like anisotropies, and lead observers to find 
that the acceleration is maximized in one direction but minimized in the opposite p^. A 
vector field may cause cosmic anisotropy and lead to a direction-dependent dark energy 
equation of state |T3]. Using Union2, Ref. [15] derived the angular covariance function of 
the standard candle magnitude fluctuations, but they did not find any angular scales where 
the covariance function deviates from in a statistically significant manner. By using 288 
SNIa [IE] , Davis et al. [Hj studied the effects of peculiar velocities, taking into consideration 
of our own peculiar motion, supernova's motion and coherent bulk motion, and they found 
that it would cause a systematic shift Aiu = 0.02 in the equation of state of dark energy if we 
neglected coherent velocities. Gupta et al. [18] introduced a statistic based on the extreme 
value theory and applied it to the gold data set of SNIa, and they showed that the data 
is consistent with isotropy and Gaussianity. In contrast to these "null- results" , Ref. [T^ 
constructed a "residual" statistic to search for the preferred direction in different slices of 
past light-cone, and they found that at low redshift {z < 0.5) an isotropic model was not 
consistent with the SNIa data even at 2-3cr. Ref. pD] found that anisotropy was permitted 
both in the geometry of the universe and in the dark energy equation of state, if we worked 
in the frame work of an anisotropic Bianchi type I cosmological model and the presence of an 



2 



anisotropic dark energy equation of state. Additionally, Refs. [2T1423] used the hemisphere 
comparison method to fit the ACDM model (and wCDM model) to the supernovae data, 
and detected a preferred axis at statistically significant level. These results are consistent 
with many other observations, such as the CMB dipole [2^, large scale alignment in the 
QSO optical polarization data [25] and large scale velocity fiows [26] . 

Ref. ^22j obtained the average direction of the preferred axes as (/, b) = (278° ±26°, 45° ± 
27°). Further analysis was made by [27], they used different low-redshift {z < 0.2) SNIa 
samples and employed the the Hubble parameter to quantify the anisotropy level, and the 
results showed that all the SNIa samples indicated an anisotropic direction at 95% confidence 
level. 

In this paper we study the plausible anisotropy in the accelerated expanding Universe 
with the Union2 data. We construct an anisotropic dark energy model and aim to detect 
the maximum anisotropy direction that deviates from the isotropic dark energy model de- 
scribed by ACDM model. Furthermore, we consider the impact of redshift on the direction 
by using the redshift tomography method, with the high-redshift Gamma-ray burst data 
as a complement. Finally, we check two other models as the description of the isotropic 
background. One is the wCDM model, the other is a dynamical dark energy model repre- 
sented by the CPL parametrization [28]. We exam that if our results are dependent on the 
isotropic background. 

The paper is organized as follows. In the next section we give a general introduction to 
the anisotropic dark energy model, which is based on the isotropic background described 
by the ACDM model, and the statistics of the model with observational data. In Sec. 
3, we give the numerical results on the maximum anisotropy directions from the SNIa data 
and high-redshift GRB data with different slices of redshift. We also give the result for the 
wCDM and CPL model. Our conclusions are presented in the last section. 

2. ANISOTROPIC DARK ENERGY MODEL 

If dark energy has anisotropic repulsive force, it will directly affect the expansion rate of 
the universe, leading to the anisotropic luminosity distance. This effect should be observable 
by the luminosity of SNIa. In this paper, we use the Union2 dataset [9J, which contains 557 
type la SNIa data covering the redshift range z = [0.015, 1.4]. 
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We try to quantify the deviation from the isotropic background as dipole modulation. 
By using the luminosity distance we define the deviation from isotropic expansion as 

dL{z) - dl{z) 



^g(z){z-n). (2.1) 



where, the true luminosity distance of the supernova is dL{z), and in an isotropic background, 
the luminosity distance is d1{z). g{z){z-n) is the modulation part of the luminosity distance, 
which makes the real luminosity function anisotropic. Note that the modulation could be 
any power-law form of {z ■ h), such as {z ■ hy, where s is a constant, but we focus on the 
dipole modulation here (s = 1). In principle, one can use this model to discuss quadrupole 
modulation (s = 2), octupole modulation [s = 4) and higher moments (s > 4). z is the unit 
direction vector of the supernova, which can be expressed by using the Galactic coordinate 
system, fi is the direction of dark energy dipole, which is the maximal expanding direction, 

n = (cos sin 6*, sin sin 6*, cos 6*) (2.2) 

where 9 G [0, tt) and G [0, 2tt). For instance, we can parameterize g{z) with linear function 
of Z as 

g{z) ^90 + giz. (2.3) 

where go and gi are two constants, representing the strength of modulation and the time 
evolution of modulation, respectively. Of course, the simplest case is where gi = 0, which 
corresponds to the case where the direction-dependent modulation is constant over all red- 
shifts. 

In a spatially flat isotropic cosmological background, the luminosity distance can be 
expressed as 

dl{z) = {l + z) f^^^dz'. (2.4) 

where Hq is the current Hubble parameter. Accordingly, the theoretical distance modulus 
^th is defined as 

Hth{z) = 5 logio d^z) +Ho , /^o = 42.384 - 5 log^o h , 

where, Hq — 100/i km ■ ■ Mpc~^. 

Since ACDM model is consistent with current astronomical observations, it is reason- 
able to take ACDM model as the isotropic cosmological background, in which the Hubble 
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parameter can be expressed as 

H\z) = Hl[n^o{l + zf + (1 - fi„o)]. (2.5) 

where is the current value of the energy density fraction of dark matter. While in the 
case of wCDM model, the equation of state of dark energy is parametrized by a constant 
w, w = p/p, therefore we have 

H\z) = H^lQmoil + + (1 - a„o)(l + ^)'+'"]. (2.6) 

And if the background is described by the CPL parametrization, the equation of state of 
dark energy is w = Wq + Wij^. Accordingly one obtains 

H{z) = H',[Q^o{l + zf + (1 - fi^o)(l + ^)3(i+'"°+'"^) exp f-^V (2.7) 



1 + 

We employ the Union2 dataset to constrain the anisotropic dark energy model. The 
directions to the SNIa we use here are given in Ref . [15] , and are described in the equatorial 
coordinates (right ascension and declination). In order to make comparisons with other 
results, we convert these coordinates to the galactic coordinates (Z,6) [22] • 

Let us suppose the experiment error between each measurement is completely indepen- 
dent, so the covariance matrix can be simplified as the diagonal component, and the xL can 
be written as 

cy^z,) 



'X.sn ^ ^ 



i=l 

where, ^"^'^{zi) is the measured distance modulus from the Union2 data, and jJ^^^Zi) is the 
direction-depedent theoretical distance modulus. 

We can eliminate the nuisance parameter /xq by expanding with respect to /xq [30] : 



xL = ^ + 2i?/xo + C/x^, (2.8) 

where 

[pth{zu /UO = 0) - ^iobs{Zi)f 



•4 = 1: 

i 
i 

c = e4 



fJ'thiZi] /io = 0) — fiobs{Zi) 



cr'^izi 



cr'iziY 



The xin has a minimum as 

xIn = A- B'/C , (2.9) 

which is independent of fiQ. This technique is equivalent to performing a uniform marginal- 
ization over /xq [30]. We will adopt xL the goodness of fitting instead of xL- 



Combining with Eq. (2.4) and substituting each anisotropic model in the xin^ can 
easily calculate the likelihood function of each parameter by performing the Markov Chain 
Monte Carlo analysis. The parameters need to be constrained are {go, gi, 9, 0), where {9, 0) 
is the direction of modulation. We then convert {9, (p) into galactic coordinate (/, b). 



3. REDSHIFT TOMOGRAPHY 

Following the procedure introduced in Sec. 2, one can obtain the best-fitting values 
and errors of parameters by performing the Monte Carlo Markov Chain analysis in the 
multidimensional parameter space. During the procedure, we use the best-fitting values 
obtained by WMAP group to calculate (fl{z) in the isotropic background. We choose 
Qmo = 0.274 for ACDM model, Qmo = 0.274, w = -1.037 for wCDM model, and 
^mo = 0.274, wq = —1.17, Wi = 0.35 for CPL parametrization, which are obtained us- 
ing WMAP + eCMB + BAO + Ho + SNe datsets 0]. 

In order to explore the possible redshift dependence of the anisotropy, we implement a 
redshift tomography analysis, taking the same procedure as before for the following redshift 
shces: 0-0.2, 0-0.4, 0-0.6, 0-0.8, 0-1.0, 0-1.2, 0-1.4. Our results for ACDM model are 
summarized in Table [Tj We also use the GRB data to explore the high-redshift behavior, 
which can extend our detection to z = 6.6. In this paper, we have compiled 67 GRB data 
in total as shown in Table [6j These samples are selected from [HI [T2] and we add in the 
position of each data point (for more details, please visit the website [25] )■ The last row of 
Table [l] shows the fitting result by combining Union2 Type la supernovae with the GRB 
samples. 

The redshift tomography analysis here shows that the preferred axes at different redshifts 
are all located in a relatively small redion of the Galactic Hemisphere (~ (/, b) = (126°, 13°) 
for Union2 data). Note that the meaning of this best-fitting direction is the same as {I, b) = 
(306°, —13°), because both directions are located on the same axis, and both are directions 
of maximum deviation from the ACDM model. The best-fitting direction is consistent with 
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redshift range 


/ [degree] 


6 [degree] 


50 


91 


- 0.2 


23-91 
^"-•^+13+212 


4-7-59 

^+30+41 


n n4t;-0.044-0.045 
'-'•"-'^^+0.026+0.063 


n 9OQ-0.713-1.648 
U-^oy+0.410+1.013 


- 0.4 


19^-18-60 
-^^^+17+255 


0-7-88 
°+24+35 


n n'^fi~0 028-0.036 
*-'■"-' "J"+0. 016+0. 025 


n nQ^-o io'i'-o.25i 

"-'•^^"5+0. 163+0.354 


- 0.6 


1 91 -18-81 
-^^^+17+255 


.Q-7-88 
'^+24+35 


n n'^'i-0 033-0.033 

U.UOO_,_Q Q]^g_^Q Q2g 


n 077-0-039-0.075 
U-'J' '+0.099+0.163 


- 0.8 


1 Qo — 22 — 88 
-'^^°+17+265 


1 c;-25-100 
-^^+19+63 


n n-^n-o oso-o.oso 

'-'■"-'"^^+0.013+0.024 


n 099-0-071-0.114 
^'■"-'^^+0.064+0.185 


- 1.0 


1 90-27-89 
i^'J+20+269 


.4-22-79 
^^+20+55 


n (Til -0.031-0.031 
U.UOl^Q Q]^2+0.024 


n n/in-0-050-o.097 

^•"J^^+0. 056+0.092 


- 1.2 


i9cr-33-90 
^^^+19+265 


1 9-29-71 
^^+20+43 


n n'?9-0 032-0.032 
"-'•"-' "J^+0. 013+0. 032 


n n/1 9-0.047-0.086 
"-'•U^^+0. 050+0.082 


- 1.4 


-.9^-26-92 
-'^^"+17+266 


1 0-25-82 
^"^+19+39 


n n'?n-o o30-o.030 
"-'•"-' "J'-'+o. 010+0. 021 


n fTil -0.039-0.075 
U.UOl_,_Q Q42_,_0.069 


- 6.6 


23-96 
-'^^^+16+293 


-1 /5-10-89 
^"+17+54 


n n9S-0-028-0.028 
"-'■"-'^°+0. 013+0. 021 


n nnR-0-034-0.055 

U.UUU_,_Q Q25_,_Q 052 



Table 1: Constraints of the directions and amplitude of maximum anisotropy for different redshift 
bins of the Union2 data and GRB data. The error-bars quoted are la and 2a errors. 



the result in Ref. [31] at la confidence level. We also re-examined the dark energy dipole 
by using the Union2 data with the same method proposed in Ref. [31], showing that the 
dark energy dipole is indeed aligned with the corresponding fine structure constant cosmic 
dipole, where the dark energy dipole direction is (/ = 309.4° ± 18.0°, b = —15.1° ± 11.5°) 
and the fine structure dipole direction is (/ = 320.5° ± 11.8°, b = —11.7° ± 7.5°). 

We plot the likelihood of the parameters {go, gi, I, b) in Fig[l] It is obvious from Eq (2.1) 
that there are two maximum anisotropic directions, which are settled on the same axis, 
accordingly g{z) can be positive or negative. Note that in our calculations, we have set 
(7o > 0, but the maximum anisotropic axis does not change if we set go <Q. 



Furthermore, we constrain the redshift dependence of the anisotropy by using SNIa data 
located only in different redshift bins: 0-0.2, 0.2-0.4, 0.4-0.6, 0.6-0.8, 0.8-1.0, thus to avoid 
influences from other bins. Accordingly, we obtain 220, 124, 102, 50 and 41 data points in 
different redshift bins respectively. Our results for AC DM model are summarized in Table |2} 
Figure [2] shows the best-fltting g{z) with la error at different redshifts. It is clear that the 
anisotropy level g{z) is close to 0, which means that the ACDM model is still consistent 
with the Union2 data very well. The error increases as redshift increases, and we can not 



7 




0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 -0.10 -0.05 0.00 0.05 

gO gl 




50 100 150 200 -40 -20 20 40 

1 b 



Fig. 1: Likelihoods for parameters {go, gi,l,b) using the fuh Union2 data, with the best-fitting 
parameters {go = 0.030, gi = -0.031, / = 126°, b = 13°). 

exclude the case g{z) = even at 68.3% confidence level. Note that the anisotropic direction 
changes as redshift bin changes, but theses directions are consistent with each other at 68.3% 
confidence level. 

By considering the case of g{z) = go and using the full Union2 data, we obtain Table |3| 
which shows that the maximum deviation from the AC DM model is 0.024, and the maximum 
anisotropy direction is similar to the results in Table [TJ but we still can not exclude the case 
(7o = at lo" confidence level. 

We also consider the wCDM and the CPL parameterized dark energy models as the 
isotropic background. We use w and (wq, wi) as the isotropic background dark energy 
parameters and fit our anisotropic parameters, respectively, and the results are summarized 
in Table |4j The results of constraints on (/, 6, go, gi) are not much different from the case 
of the ACDM model. This means that the best-fitting value of the maximum deviation 
direction from the isotropic background is not sensitive to the details of isotropic dark 
energy models. 
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redshift range 
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00/1-279-301 
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-^°+53+108 


n /1 87-0.487-0.487 
+0.412+0.511 


n 077-0.039- 0.075 
"-'•^' '+0.099+0.163 


0.6-0.8 


lzlQ-115-116 
-•■^^+241+244 


f;q-106-140 
""^+10+25 


n icq-0.163-0.163 
U.lUO_^Q 78y_,_0 


n 1Q1 -1.138-1.257 
"-'•-^^-'-+0.270+0.330 


0.8- 1.0 


96-113 
-■-^"+233+246 


/.-33-68 
"+46+90 


n 007-0.996-0.997 
'-'•^-^'+0.003+0.003 


1 -11-1-0.293-0.293 
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Table 2: Constraints on direction and strength of modulation for several redshift bins of the Union2 
data fitting with the ACDM model, together with the la and 2a errors of parameters (/, b, go, gi). 

1 ' ' : 

0.2 0.4 0.6 0.8 1.0 1.2 

z 

Best-fitting g{z) with la error at different redshift. 



/[degree] 


6 [degree] 


50 


-190-23-94 

-'-^"+22+264 


14_97 
-•-"+21+64 


n 09/1 -0 024-0. 024 
"-'•"-'^^+0.008+0.021 



Table 3: Constraints on direction and strength of modulation for the full Union2 data fitting with 
the ACDM model, together with the la and 2a errors of parameters (/, b, go). 

In the end, we make Bayesian statistical comparison between the isotropic dark energy 
models and the anisotropic dark energy models, in order to see which model is more favored 
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Fig. 2: 
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Models 


I [degree] 


5 [degree] 


50 


91 


wCDM 


-1 OO — 28 — 88 
-'-^°+17+256 


30-74 
-^^+17+39 


n no/i -0.033-0.034 
"-'•^"-'^+0.013+0.026 


n n/ir^-O-056-0.093 
^•^U+0.049+0.090 


CPL 


1 97-21-80 
+17+230 


1 c-32-77 
-^^+17+61 


n n'^f;-0-035-0-035 
"-'•^"-'"-'+0.015+0.027 


-n 040-0.046-0.347 
"-"•"-'^^+0.049+0.102 



Table 4: Directions of maximum anisotropy fitting with the wCDM model and CPL model, 
together with the la and 2a errors of parameters (/, 6, 9i)- 

by observational data. The Bayesian Evidence E provides us with a good metric to quantify 
the level of consistency between each model and observational data [53] . 

BE = j C{D\9,H)p{e\H)d9. 

The ratio of Evidences between the two models Bab = E{A)/ E{B), also known as the 
Bayes factor, measures the goodness of fit of the models. In the following, we calculate 
the Bayesian Evidence ratio of the anisotropic model (A) and the isotropic model (B) by 
assuming the uniform prior for every parameter. According to the Jeffreys grades, the result 
is summarized as follows (see Table [s]). 



Models 


Parameters 


Bayes factor 


Comparison 


ACDM 


^mO 


-0.35 


Weak 


wCDM 


f^mO, W 


1.08 


Significant 


CPL 


^mOi Wo, Wi 


3.17 


Strong to very strong 



Table 5: Comparison result using the Bayes factor for anisotropic and isotropic dark anergy models. 



(1) If the background is described by the ACDM model, the parameters for isotropic 
and anisotropic models are (fimo) and (/, b, go, gi), respectively. The Bayes factor is 
—0.35, which means that there is "weak" evidence that the ACDM model is better 
than the anisotropic model. 

(2) If the background is described by the wCDM model, the parameters for isotropic and 
anisotropic models are {Qmo, w) and (Z, b, go, gi), respectively. The Bayes factor is 
1.08, which means that there is "significant" evidence that the anisotropic model is 
better than the wCDM model. 
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(3) If the background is described by the CPL parametrization, the parameters for 
isotropic and anisotropic models are {VLmQ, Wq, Wi) and (Z, 6, Qq, gi), respectively. The 
Bayes factor is 3.17, which means that there is "strong to very strong" evidence that 
the anisotropic dark energy model is better than the CPL parametrization. 

A. Distinguish Anisotropic Dark Energy With Peculiar Velocity Field Model 

In this paper, we consider the model of anisotropic dark energy, which may potentially 
be degenerated with the peculiar velocity field model in an isotropic background. For ex- 
ample, the best-fitting bulk flow direction found in Ref. [M] is quite close to the direction 
(/, b) = (306°, —13°). . This is because, the peculiar motion of galaxies or supernovae can 
also produce a direction-dependent luminosity distance, which may indicate some preferred 
direction. With the current data at low redshift, we are not able to distinguish them, but 
we propose the following method which is useful for doing this in advance of new data. 

First, if the anisotropy is caused by the peculiar velocity, the anisotropic direction should 
be randomly distributed on different cosmic scales, because peculiar velocity is driven by 
emergent of large scale structure, but if the anisotropy is caused by the dark energy dipole, 
the anisotropic direction should be a constant on all cosmic scale, due to the non-local effect 
of dark energy. So redshift tomography method may tell the differences between the two 
models if high-z SNIa data are available. 

Second, peculiar velocity is a local effect and should be zero if averaged on the whole 
cosmic scale [32]. But the dark energy dipole should not change with the redshift. So by 
average the galaxy luminosity over a large volume, one can distinguish where the direction- 
dependence of luminosity is due to peculiar velocity or dark energy dipole. 

In addition, the Integrated Sachs- Wolfe effect (ISW) can be used to detect the dark 
energy dipole, because if the accelerating expansion is anisotropic, the photons that travel 
from different distances in different directions are potentially observable. On the other hand, 
peculiar motion of galaxies can only re-scatter the CMB photons and produces secondary 
anisotropic effect, which acquires its maximum at galaxy and cluster scales. 
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4. CONCLUSIONS 



There is a tentative evidence that the anisotropic direction on the cosmic expansion 
exists. If such a cosmological preferred axis indeed exists, one has to consider an anisotropic 
expanding Universe, instead of the isotropic cosmological model. 

In this paper, we investigated the plausible anisotropy in the accelerating expansion 
Universe with the Union2 data and high-z GRB data. We construct an anisotropic dark 
energy model, where we quantify the strength of modulation as function g{z) — go + giz and 
direction of maximum expansion as {I, b). We found that the maximum anisotropic deviation 
direction is {l,b) — (126°, 13°) (or equivalently (Z, 6) = (306°, —13°)), and the maximum 
anisotropy level is described by the parameters go — 0.030^o!oio> fl'i — ~0-031+o!o42 (obtained 
using Union2 data, at la confidence level). Furthermore, we used the redshift tomography 
method and found that the anisotropic direction depends on the redshift. However, it is clear 
that the anisotropy strength g{z) is close to within la confidence level, which indicates 
that there is no strong evidence against isotropic AC DM model. 

If the modulation is constant over all redshift, by using the full Union2 data, we found 
that the maximum deviation from the AC DM universe is 0.024;^:°^^, and the maximum 
anisotropy direction is similar to the former results, but we still can not exclude the case 
go = aX la confidence level. 

We also discussed the cases where the dark energy equation of state is described by a 
constant w and w{z) — Wo + ^i^^^ respectively, but the results show a similar anisotropic 
direction. This indicates that the best-fitting value of the maximum deviation direction 
from the isotropic background is not sensitive to the isotropic dark energy models. 

Finally, by using the Bayesian Evidence, we found that the anisotropic dark energy model 
docs not show great statistical evidence better than isotropic wCDM model, except that 
there is a slightly greater evidence of anisotropic dark energy than GPL isotropic dark energy 
model. 
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Table 6: Selected 67 GRB data points, added in positions described in the equatorial coordinates 
(right ascension and declination). 
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